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A STEERING LAW FOR THREE DOUBLE-GIMBAL CONTROL 
MOMENT GYRO SYSTEMS 

I. INTRODUCTION 


Control moment gyro (CMG) systems are effective actuators for attitude 
control of a large variety of spacecraft. In Skylab, a doublc-gimbai CMG system 
was used successfully. Also, in some future projects, like the Large Space 
Telescope (LST) and the High Energy Astronomy Observatory (HEAO), CMG 
systems are proposed as the primary system for attitude control. 

In this report, a new steering law, whose principle is very simple, for 
three double-gimbal CMG systems is proposed. It is applicable to almost any 
configuration, and the CMG-out operation requires no special consideration. 
Although several steering laws for double-gimbal CMG systems have been 
developed in the past including one that was used in the Skylab [ l-Sl , the basic 
principles of these steering laws are more complicated than that of the proposed 
steering law. 

The main idea is utilization of a formula for the general solution to a set 
of linear algebraic equations whose number is less than the unknown variables. 
The general solution is given in terms of the pseudoinverse of the coefficient 
matrix for the set of linear algebraic equations and it contains an arbitrary 
vector. This arbitrary vector is used to obtain a desirable momentum distri- 
bution among the CMGs. 

In Section II, a general description of the steering law is given. In 
Section III, two examples are worked out. The first example is a three double- 
gimbal CMG system with an orthogonal configuration and the second one is a 
three double-gimbal CMG system with a parallel configuration. In both 
examples, digital simulation results are given to illustrate the effectiveness of 
the steering law. 



II. DESCRIPTION OF THE STEERING LAW 


The output torque of any CMG system can be described by 


6 = C(M 6 


where 5 is the n-dimensiot.^1 gimbal angle vector (in the case of three double- 
gimbal CMGs, n= 6), £ is tl . three-dimensional output torque vector, h is 
the three-dimensional total angular momentum vector, and C(5) is a coeffi- 
cient matrix, each entry of which is a function of 5 . h and b denote the 
time derivatives of h and ^ , respectively. It is assumed there are no 
gimbal stops on the CMGs. 

Given the present state 6 >£ the system, it is required that one obtain 
a gimbal rate command r which satisfies 




( 1 ) 


where is the command torque. The general solution to the above equation 
is given by [ 4 ] (C(6) is also denoted by C hereafter] 


r 


C*L + II-C*C]k 
c — 


( 2 ) 


where C is the pseudoinverse of C , I is the n x n identity matrix and k 
is an arbitrary constant 'sector. The second term can be used for momentum 
distribution . 

It is assumed that a proper scalar criterion function which measures the 
desirability of the momentum distribution of the present state 5 is given, 


w = f(6) 


(3) 


2 



and it is also assumed that the larger the value of w , the better the momenb<m 
distribution. The time derivative of w is given by 


w = 4 0 


where the superscript T denotes transpose and 


4 4 1 

— 12 n 


f = iM) 




(4) 


Under the as j”mption that the gimbal rate command ( 2) is perfectly realized, 
that is, r, the time derivative of w is given as follows: 


w = 


4^C^t + 4^tI-C*Clk 

"C — — 


(5) 


In order to increase the value of w , it will be best to select k as 

k = 4kj (6) 

where k^ is a positive constant. The reason for the selection of equation (6) 
follows: From equations (2), (s) and (6), one obtains 

r = C*t + [I-C*Cl i k, (7) 

w = + 4^lI-C*Cl 4kj . (8) 

and the second term in the right-hand side of equation (s) becomes nonnegative 
§ 

(because (I-C C] is idc'^potent) , contributing for an increase of the value w. 


3 
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Equation (7) is the basic equation of the steering law. It involves the 
selection of only one scalar constant . 

In summary, the main procedure of the steering law is given as follows; 

1. Calculate C . 

2. Calculate 4 • 

3. Calculate r by equation (7) . 

Remark 1 ; The meaning of the constant k^ will be considered. If 
t ^ = 0 ( zero vector) , then 


r = r I = (I-C C] |_k 

t = 0 

"C — 


( 9 ) 


w = w 1 = ^^[I-c’^Cl 4_k 

” 't = 0 ^ 

-c — 


( 10 ) 


The maximum of w under the constraint Hr II :£ 1 ( II II denotes norm) is 
given by 


w 

Omax 




( 11 ) 


Hence from equations (9) and (ll) 



1 Omax 


( 12 ) 


Therefore, if = £» the magnitude of r becomes proportional to and 

kj is the proportional coefficient. In practical applications, in order to 


4 


prevent from becoming too large, a limitation for should be given. 

One candidate for this is Hr ^11 s k„r , i.e. , 

-0 2 g 

\ ^ k^r^ (13) 

where r (deg/sec) is the hardware limit for tlie gimbal rate and k is a 
S 2 

constant (O k^ £ l). 

T 

Remark 2; If det(CC ) =^0, then 

C = C (CC ) . (14) 

On the other hand if det(CC ) = 0, the calculation of C is rather compll- 
cated. A simple way to overcome this difficulty will be to approximate C by 
Cp where Cp=C(^+^) and is a small perturbation of 6 juch that 

det(C Cp) 4 0. For double-gimbal CMC systems dct(CC ) = 0 happens 

only when all the momentum vectors of the CMGs are on the same line (parallel 
or antiparallel) . Such a situation can hardly happen except when n = 4 (i. e. , 
two double-gimbal CMGs). Also one of the main purposes of the momentum 
distribution is to avoid such a situation. Hence the small error caused by the 
# # 

approximation of C by Cp will be negligible in practical applications. 

T 

Remark 3; In the following section, the root of the determinant of CC 
will be taken as the criterion fimction f( 6 ) ; 

f(6) = n/ detG 


where 


G = CC 


Igijl 


l,j = 1, 2, 3 


5 



~S3*1 . . 


I 


The reason for this selection of f(6) will be given later. For this criterion 
fucctiont is given by 


t = ^ y g. d f = 1 2 f 

® ®i 2-«ra5G i,j«i ' '* ‘J 


( 15 ) 


where 


I'ij ' ’ 


‘ Vkk - «ik 


(16a) 


~ ~ ^ij^kk * ~ (lf2f3)» (2|3(l) or (3yl|2) . 

(16b) 


Using the fact that g = g.. and d d . equation ( 15) can be rewritten as 
follows: ^ ^ 


h **11 l®22 **22 i®33 **33^ * 1^12 ^12 ^ 1^*13 **13 


+ ,g’ d J/*^detG 
1*^23 23 


(17) 


III. EXAf^PLES 

In this section, two examples of three double-gimbal CMG systems will 
be given to show the effectiveness of the steering law. 

Example I t For a three double>gimbal CMG system in an orthogonal 
configuration, several notations will be introduced. For the ith CMG (i = 
1,2,3), an orthogonal coordinate system (X^ Z^) is used, where axis is 


I 



the outer gimbal pivot axis and and axes are fixed with respect to the 

spacecraft, a, and are the outer and inner gimbal angles, respectively, 

(in degrees) of tiie ith CMG, as shown in Figure 1. The angular momentum 

3 

vector of the ith CMG is denoted by h and the total momentum by h(- ^ h .) . 

^ ” i=l ‘ 

h^ denotes the magnitude of h.. Then, inX^Y.Z^ coordinate system. 


h, ca. 
ill 


hi= \ 


h s0 
i^i 



OUTER OIMtAL PIVOT 


Figure 1. Outer and inner gimbal angles of a double-gimbal CMG. 
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and 


h. = 
—1 


-h. Ba, c/3. - h. cor, sfi, 

ill ill 




h, Q^. 

i 1 


or. 


L^i 


where c6 = cos6 , s5 = sin6 . 

Now the three double<^imbal CMG system in an orthogonal configuration 
shown in Figure 2 is considered. Let 


5 


2i-l 


a. 


1,2,3 


621 = i = ^* 2,3 


then, in XYZ -coordinate system given in Figure 2, 


C = 


-h sor c0 -h cot bB 


h cB 
2 2 


hcoc^ -hscrs3 -hsorc^ -hcos/3 

2 2 2 2 2 2 


“i “^1 \ “2 “^2 -“a ““2 "^2 


"3 “3 "^3 -"3 ”3 ”^3 

» ■'3'^3 

-^3 “3 “^3 -^3 “3 ^3 
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Figure 2. A three-CMG system in an orthogonal configuration. 


Let G = CC^ and f(0) = »>/detG", then ^g, 

2Ul*ll “ ‘2Ul®jj " 

2Ul®iJ “ ’ 2 uA ' 

21- A ' *1 ’ 2l-l*kk 


*8 are given by 


*5 “i % 

= 0 



sc' * - k' = 2tf b a 
2^11 2 ni I al i 




SC* = 4tf a a_ . 2' = - tf ca 

2rij I al * 2rik i i fil 


2i*^Jk “ ’ 2i*kk " “ 


^l 


(18b) 


where (i»j,k) « (1,2,3), (2,3,1) or (3,1,2), and 


a ~saca . - 80 c0 

al i i * 01 


b = <^a - 8*a , b- = c*8 - s*0 

ai i i ’ 0i 


(19) 


Therefore, 4^*s are given by equations (16), (19), and 


4 = h*{[a (d - d ) - b d I b 

^21-1 r al' ii al Ij j 


01 


- Isa, d„ + ca, d„ I a.,}/'v/detG 
i ik i jk 01 


(20a) 


4 “ h*( [b (d -d)-d +4a d]a* 

^2i r ar ii jj' kk a I Ij ^l 

' I ^k ■ “'i V 


(20b) 


Hence the steering law can be described as follows: 
1.. For a given 6, calculate C(5) and detG. 

2. If detG » 0, then let 


10 


i 



6^ = 6^ + A6 , I = 1, 2, .... 6 


anc go to step 1. K detG ^ 0 » go to step 3. A6 is a small positive number 
(typically A6 » 0.2). 


^ T —1 
3. Calculate C = C G” . 


4. Calculate from equations (16), (19), and (20). 


5. Calculate r^=C t . 

— t “C 


If the absolute value c£ any component of r 


is larger than r , limit each conq>onent of r proportionally, such that die 
maximum absolute value of conqionents of r ^ is equal to r^ . 


6 . 

where 


k 


1 


Calculate r = r + [I-C C] ^ k 
= min (kg , k^r^ {£^[I-C*Cj g ) 


and k and k are positive constants. If the absolute value of any component 

A u 

of r is larger than r , k, is decreased until the absolute value of every 
- g 1 

component of r becomes smaller than r^ . 


This completes one cycle: After one unit time interval A, the state of 
the system 5 is measured and the next cycle begins. In normal operation, the 
qualifying remarks mentioned in steps 5 and 6 will very seldom occur. 


There are eight zero momentum stationary states. One is given by 
{o!^ = -45, /3^ = 0 , i = 1 , 2 , 3 } as shown in Figure 3, and die state which is 

symmetric to this state with respect to the origin is also a zero momentum 
stationary state. Others are symmetric to these two states with respect to the 
XY -plane, the YZ -plane, and the ZX-plane. In Figure 3, the outer circle shows 
the contour of the sphere with a unit radius and three inner circles denote the 
cross sections of the sphere by the XY-, YZ - and ZX-plane. The three arrows 
terminating on the sphere surface denote the three momentum vectors h 
i» 1,2,3. 
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I 

I 



Figure 3. Zero momentum stationary state of the three-CMG 
system in an orthogonal configuration. 

A digital simulation has been completed using a mini -computer. Values 
of constants were selected as follows: 


h^ - 1 , i - 1,2,3, (normalized) 

— 0. 2 , kg = 0. 1 

r = 2,A=8,A5 = 0.2 
g 




} 



Figure 4 shows a response of the system to a torque command in the X axis 
direction with magnitude 0.01. The initial state is the zero momentum state 
shown in Figure 3. Dots on the unit sphere show the end points of the momen- 
tum vectors h. , i = 1»2»3, and arrows along the side of the dots show the 

time trajectory (broken lines indicate that the vectors are located on the back 
half of the sphere). The output torque followed the torque command 
satisfactorily. 



Figure 4. Response to an X directional torque command with 
magnitude 0. 01 ( orthogonal configuration) . 

If the direction of the command is exactly the opposite of a vector h^ 

all through the period of steering, there is a possibility for the system to get 
intr. a singular state. For example, if the initial state is = -45 , = 0 , 


i 

\ 


i= 1,2»3} and the torque command is (t = 0.005, t = 0, t = -0.005} (see 

X ^ z 

Figure 5) , then h does not move at all and the state falls into a singularity. 


In practice, however, such a situation can hardly happen because of the exis- 
tence of disturbances on the spacecraft and noise in the sensors. Figure 6 
shows a response to a torque command in the neighborhood of the most um'avor- 

able direction. The initial state is {a = = -45, a - -44, /3 = 0, i = 1,2,3) 

X 2 o 1 


auu the torque command is {t = 0.005, t = 0, t = -0.005}. It is seen from 

X y z 

the fi^re that by the compensation of h and h , h is turned to a desirable 
direction. 



Figure 5. Possible direction of torque command for the system to 
get into a singular state ( orthogonal configuration) . 




i.i7v 
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Figure 6. Momentum distribution process for an unfavorable torque 
command (orthogonal configuration). 

Moreover, even if the system falls into a singular state, the steering 
law has an ability to recover from it, as shown in Figure 7 where the initial 
state is {a = ^ = a = /3 = 0 , /3 - '” = -90) and the torque command is kept 

11^0 £t ij 

at zero. For a C MG -out operation, the only change that should be made is to 
set h^ = 0 when the ith CMG has failed. The control performance of this 

system was very satisfactory in every case simulated. 


Example 2; A three double-gimbal CMG system in a parallel configura- 
tion. The same notation used in Example 1 is valid in this example. The three- 
CMG system in a parallel configuration* as shown in Figure 8 is considered. 

For this configuration, the matrix C(6) is given by 


1. Kennel, H. F. : Control Law for Parallel Mounted Double-Gimbal Control 
Moment Gyros. Second Annual Research and Technilogy Review, Marshall 
Space Flight Center, Oct. 1974. 


15 




b. Positions of momentum vectors in XY-plane at time r . 
Figure 7. Recovery from a singular state (orthogonal configuration). 



1 



I 


s^’T- 



Figure 8. A three-CMG system in a parallel configuration. 


-hjSajOUj -hjCKjafij ->> 2 <» 2®^2 

C(6) - I hjCOfjCSj I'2'“2'^2 -*2 ”2*^2 


hjCSi 


"2<^2 


-h sa„ c/3„ -h c0!„ s^„ 
3 3 3 3 3 3 


h„ CQ' c^„ -h sof- s^„ 
3 3 3 3 3 3 


h„ c/J„ 

3 3 
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Let £(6) - N/detG , then by a derivation similar to that in Example 1, i is 
is given by equations ( 16) , ( 19) , and ^ 


1 = h^{ta . (d_ -d ) - b d ] b 
2i-l 1*^ ai 11 22' Ofi 12 0i 






01 


- '“l ‘*13 - ®“l V 


(21a) 


(21b) 


The procedure of the steering law is the same as that in Example 1 with 
the following exception; Use equation (2l) instead of equation (20) in step 4, 
and instead of step 2 the following step should be used to ensure the avoidance 
of the case of detG = 0 : 

If detG = 0, then let 


a; = a, + A 6 *(i .3.5) 


^ = 1, 2, 6 


and go to step 1. If detG #0 go to step 3. Aa* is a small positive number. 
Zero momentum stationary state is given by (la -a| = |a -a I 

II 1 12 2 3 

= la.^ -a^l -- 120, 0.= 0, i= 1,2,3). One such state is shown in Figure 9; 

any other state can be obtained by a rotation of this state around the Z axis. 
A digital simulation was performed with 


= 0.2 , kj = 0.1 


= 2 , A = 8 , Aa' = 0.1 




IS 

iVALlTY 
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Figure 9. Zero momentum stationary state of the three-CMG 
system in a parallel configuration. 

Figure 10 shows a typical behavior of the system; a torque command = ly = 

t = 0. OOS) is applied to the system with a zero momentum stationary state 
z 

{o' = 0, Q! = 120, Of = -120, /3. = 0, i- 1,2,3). Figure 11 shows a response 

id O 1 

to a torque command in the vicinity of the most unfavorable direction. The 
initial state is selected as (a = 60, a = 179, a = -60, /3. = 0, i= 1,2,3), 

and the X-directional torque command with magnitude 0*005 is applied* The 
figure shows a very good momentum distribution procedure* 

If the initial state is {oi = 60, a = 180, a ~ -60, = 0, i -- 1,2,3} 

} d O 1 

and the torque command is kept exactly in X axis direction, the system will f 'll 




Figure 10. Response to torque command t = t = t = 0.005 

X y z 

(parallel configuration). 

into the singular state shown in Figure 12(b). r - 0 . As mentioned in Example 

1, however, the possibility of occurrence of this situation is almost zero. 

Moreover, even if such a situation does occur, the steering law has the ability 

to recover from it. as shown in Figure 12 where the initial state is {Oj = “ 0* 

a = 180, 0 - 0. i = 1 . 2 , 3 ) and the torque command is kept at zero. 

3 1 


IV. DISCUSSION 

Although the root of the determinant of CC (-*^detG) has been used as 
the criterion function, any reasonable funct ion of 6 can be a candidate for 
it. The reason for the selection of *JdetG in this report is as follows: 
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1 





b. Positions of momentum vectors in XY -plane at time r. 

Figure 11. Response to an X directional torque command 
(parallel configuration). 
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■ 0 r - 264 

b. Positions of momentum vectors in XY-plane at time r . 

Figure 12. Recovery from a singular state (parallel configuration). 
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1* detG is a scalar variable which is most directly related to the 
singularity of the system. 

2. By digital simulation, it was foimd that VdetG gives £ which is 
large enough to drive the system out of the vicinity at the singular states, while 
detG gives very small around the singular states. 

3. So far, the author co uld no t find any simpler criterion function which 
gives as good performance as *^detG . 

For some specific CMG configurations, there may be simpler criterion func- 
tions. For example, for the parallel conffgur.tion given in Example 2, 




is a good candidate for f(6), whp"' o;^^, i= 1,2,3 and are reference 

angles that should be specified by the designer as functions of 6 . Whether or 
not this is really good for f(5) could be determined only after investigating 
various methods of giving and and performing a digital simulation 

study. To find such a simpler criterion function is a topic for future 
investigation. 

•*% 

The author believes that the software of the proposed steering law is 
much simpler than any other steering laws developed so far. Rigorous com- 
parison study should be performed, however, to give any objective statement 
about computer software. This is another direction of the future study. 

Since ai^ r which satisfies equation (l) can be expressed in the form 
of equation (2) by selecting an appropriate k, ai^ steering law which produces 
the same torque output as the torque command could be expressed in the form 
of equation (2). Also, if k is taken to be zero, equation (2) becomes a 
steering law which has usually been called the pseudoinverse steering law. 
Therefore, the steering law developed here may be called the "universal 
pseudoinverse steering law." 

Systems utilizing only three double-gimbal CMQs have been developed 
in some detail. However, it is clear from the general discription given in 
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Section 2 that the steering law has the possibility of remaining elective for any 
number of double-gimbal and/or single^imbal CMGs with any configuration 
(which result in n ^ 3), possibly with a proper selection of the criterion func- 
tion. Reaction wheel systems are also expressed by ( l) with a constant 
coefficient matrix C instead of a state -dependent matrix C(^). Therefore, 
the approach taken in this report will be useful in obtaining steering laws for 
reaction wheel systems. Research in this direction should also be made in the 
future. 


V. CONCLUSION 


A steering law, whose principle is very simple, for three double-gimbal 
CMG systems is proposed. This steering law is applicable to systems with 
alr.'.ost any configuration of CMGs and the CMG -out operation needs no special 
modification. 

To illustrate the effectiveness of the steering law, two examples of three 
double-gimbal CMG systems in an orth<^onal conf^uration and in a parallel 
configuration were worked out in detail. Simulation results have shown that 
any command torque can always be met except when the system is in a singular 
state and that, whenever the system is in, or close to, a singularity, the 
steering law drives the system out of the vicinity of the singularity. 
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